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Abstract 

We derive estimates for the solutions to differential equations driven 
by a Holder continuous function of order /3 > 1/2. As an application we 
deduce the existence of moments for the solutions to stochastic partial 
differential equations driven by a fractional Brownian motion with Hurst 
parameter H > i. 



1 Introduction 

We are interested in the solutions of differential equations on R m of the form 

x t =x + / f(x r )dy r , (1.1) 
Jo 

where the driving force y : [0, oo) — > R m is a Holder continuous function of 
order f3 > 1/2. If the function / : R d — > R md has bounded partial derivatives 
which are Holder continuous of order A > 4 — 1, then there is a unique solution 
x : W n — > R, which has bounded ^-variation on any finite interval. These 
results have been proved by Lyons in j^j using the p-variation norm and the 
technique introduced by Young in [H] . The integral appearing in i|l.l[) is then a 
Riemann-Stieltjes integral. 

In [7] Zahle has introduced a generalized Stieltjes integral using the tech- 
niques of fractional calculus. This integral is expressed in terms of fractional 
derivative operators and it coincides with the Riemann-Stieltjes integral fdg, 
when the functions / and g are Holder continuous of orders A and /x, respectively 
and X + fi > 1 (see Proposition ^ . II below) . Using this formula for the Riemann- 
Stieltjes integral, Nualart and Ra§canu have obtained in (2) the existence of a 
unique solution for a class of general differential equations that includes i|l.l|) . 
Also they have proved that the solution of (|1.1|) is bounded on a finite interval 
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[0, T] by Ci exp(C2 ||y|lo,T,^)' where n > ^ if / is bounded and n > 1 _ 1 2/3 is / 
has linear growth. Here ||y||o,r,/9 denotes the /^-Holder norm of y on the time 
interval [0, T]. These estimates are based on a suitable application of Gronwall's 
lemma. It turns out that the estimate in the linear growth case is unsatisfactory 
because k tends to infinity as f3 tends to 1/2. 

The main purpose of this paper is to obtain better estimates for the solution 
Xt in the case where / is bounded or has linear growth using a direct approach 
based on formula l|2.8|l . In the case where / is bounded we estimate sup 0<t<T \xt | 

by 

c(i + ll< r)/3 ) 

and if / has linear growth we obtain the estimate 

Ciexp(c 2 ||y||| )T>/3 ). 

In Theorem 13.11 we provide explicit dependence on / and T for the constants C, 
C\ and C2. 

Another novelty of this paper is that we establish the explicit dependence 
of the solution Xt to i|l.lfl on the initial condition xq, the driving control y 
and the coefficient / (Theorem 3.2). Similar results are obtained for the case 
l/3</3<l/2ina forthcoming paper 0. 

As an application we deduce the existence of moments for the solutions to 
stochastic partial differential equations driven by a fractional Brownian motion 
with Hurst parameter H > | . We also discuss the regularity of the solution in 
the sense of Malliavin Calculus, improving the results of Nualart and Saussereau 
H], and we apply the techniques of the Malliavin calculus to establish the exis- 
tence of densities under suitable non-degeneracy conditions. 

2 Fractional integrals and derivatives 

Let a, b <E R with a < b. Let / E L 1 (a, b) and a > 0. The left-sided and 
right-sided fractional Riemann-Liouville integrals of / of order a are defined for 
almost all x £ (a, b) by 

I? l+ f{t) = -±- f \t-sT~ 1 f(s)ds 

and 

ISLf(t) = t^- [\s-tr- 1 f(s)ds, 
r(a) J t 

respectively, where (—1) a = er na and T (a) — / °° r a ~ 1 e~ r dr is the Euler 
gamma function. Let I" + (L P ) (resp. I^_{L P )) be the image of L p (a,b) by the 
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operator I% + (resp. If / G I£ + (£ p ) (resp. / £ /£_ and < a < 1 

then the Weyl derivatives are defined as 

Q+ W r(l-a) \^(i-a) a y o y 1 7 

and 

Db - f{t)= w^){w^r +a L ( S -tr^ ds ) (2 - 2) 

where a < t < b (the convergence of the integrals at the singularity s = t 
holds point-wise for almost all t G (a, b) if p = 1 and moreover in L p -sense if 
1 < p < oo). 

For any A G (0, 1), we denote by C x (a,b) the space of A- Holder continuous 
functions on the interval [a, b] . We will make use of the notation 

II II \X T XQ\ 

■ a<e<r<6 r _ y r 

and 

IMU,f>,°° = sup |a; r |, 

a < r < b 

where x : R d — ► R is a given continuous function. 
Recall from [5] that we have: 

• If a < - and q — i— - — then 



=igL(L*)c£« (a.b) 



• If a > - then 
p 



I« + (LP)U IZ_(LncC a ~Ha,b). 

The following inversion formulas hold: 

Ia+{D a a+ f) = /, V/ G I " + (£*) (2.3) 

Ia-(D^_f) = /, V/eCP) (2.4) 

and 

DZ_(l%_f)=f, VfeL'(a,b). (2.5) 
On the other hand, for any /, g G L l (a, b) we have 

" I c : + f{t)g{t)dt = (-!)<* f(t)I?_g(t)dt , (2.6) 



and for / G (L p ) and g G I%_ (LP) we have 

D a a+ f(t)g(t)dt = (-I)"" / f(t)DZ_g(t)dt. (2.7) 
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Suppose that / £ C x (a,b) and g G C^(a,b) with A + /i > 1. Then, from 
the classical paper by Young [H|, the Riemann-Stieltjes integral J fdg exists. 
The following proposition can be regarded as a fractional integration by parts 
formula, and provides an explicit expression for the integral J fdg in terms of 
fractional derivatives (see [7]). 

Proposition 2.1 Suppose that f G C x (a,b) and g G C^(a,b) with A + \i > 1. 
Let A > a and /i > 1 — a. Then the Riemann Stieltjes integral J fdg exists and 
it can be expressed as 

f fdg = {-l) a f £>«+/ (t) D\z a g b - (t) dt, (2.8) 

J a J a 

where g b - (t) = g (t) - g(b). 



3 Estimates for the solutions of differential equa- 
tions 

Suppose that y : [0, oo) — > R™ 1 is a Holder continuous function of order f3 > 1/2. 
Fix an initial condition xq G M. d and consider the following differential equation 

x t = x + / f(x r )dy r , (3.1) 
Jo 

where / : R d — > W md is given function. Lyons has proved in |2] that Equa- 
tion (|3.1|l has a unique solution if / is continuously differentiable and it has a 
derivative /' which is bounded and locally Holder continuous of order A > -g — 1. 

Our aim is to obtain estimates on x t which are better than those given by 
Nualart and Ra§canu in 

Theorem 3.1 Let f be a continuously differentiable such that /' is bounded 
and locally Holder continuous of order A > 4 — 1 . 

(i) Assume that f is also bounded. Then, there is a constant k, which depends 

only on (3, such that for all T , 

sup \x t \<\xo\ + kT\\f\\ 00 \\f'\\tf\\y\\t T p- (3.2) 

0<t<T ' 

(ii) Assume that f satisfies the linear growth condition 

\f(x)\ <a + ai \x\, (3.3) 

where ao > and a\ > 0. Then there is a constant k depending only on 
[3, such that for all T, 

sup \x t \ < 2 kT ^'^ a ^ /fi ^o^(\xo\ + 1) . (3.4) 

0<t<T 
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Proof. Without loss of generality we assume that d = m = 1. Assume first 
that / is bounded. Set \\y\\p — \\y\\o,T,p- Let a > 1/2 such that a > 1 — (3. 
First we use the fractional integration by parts formula given in Proposition ^. II 
to obtain for all s,t G [0, T], 



f(x r )dy r \< J \D° + f(x r )D^ a y t _(r)\dr. 

From H2.2fl and l|2.1|l it is easy to see 

\D]Z a y t -{r)\ < k\\y\Up\t ~ r\ a+ ^ < k\\y\\ p \t r^" 1 (3.5) 

and 

\D? + f(x r )\ < k [H/IU (r - 8)~ a + H/'IL \\x\UtAr ~ sf~ a ] . (3.6) 
Therefore 

f(x r )dy r \ < kWyyjjWfW^ir-sy^t-rr^- 1 

+ ll/'ILII»II.A/j(r-«)' J - (*-r) a+ ' s - 1 ] dr 

< *lli/ll/j[ll/IL(*-^ + ll/ , ILW-,t,/3(*-«) 3/s ] • 

Consequently, we have 

lklU,^<fc||y|l4ll/lloc + ll/ / llocll^lU,t,/3(«-^] • 

Hence, 

\\x\\ s , t ,e<k\\y\\ H/ll^ (l-fcH/'ILHyll^t-^)- 1 . 

Therefore, 

IMUt,oo < |ar a | + ||x|| s , t „g(i - s) 73 

< W + ll/IL (1 ~ ll/'lloo ~ (* - »?■ 

Let A :— k Wf'W^ \\y\\p- Divide the interval [0, T] into n — T/A subintervals and 
apply the above inequality on the interval [0, A], [A, 2 A] and so on recursively 
to obtain 

sup \x t \ < \x \+kT 11/11^ \\y\\ p (l -A^y 1 ^- 1 . 

0<t<T 



With the choice A = (j-jr'j ^ we get 



sup \x t \ < lxol + ATH/ILIIyll/3 1 -^r(k H/'IL hh)^ 

o<t<r fi{l-(3) — 



= \M + kT\\f\\^-\\y\\}. 
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This proves the inequality 13. 2 [I . 

Assume now that / satisfies (|3.3[1 . In this case, instead of l|3.6(l we have 

\D« + f(x r )\ < k [(oo + ai \ Xr \) (r - S )- Q + 1 1 /' 1 1 oo INU^ " ^1 ■ 

As a consequence, 

IMkt.fl < k\\y\\ p [o + ai ||i||, )ti00 + H/'IU ||i||. ) t, / j(t - s)' 3 . 

Or 

\\x\\ s ,t, < k\\ y y (a + 01 |k|| M)00 ) (1 - k ll/'l^ - sfy 1 . 
Therefore, 

n < ix^+kwyyii-kwrw^wvut-sfy 1 

x (a + ai INI,,*,*,) ft-*)"- 

As before, divide the interval [0, T] into n — T/A subintervals and set A = t — s. 
Denote 

A = k\\f 



We have 



This implies 



ioo \\y\\p 
B = ka \\y\\p 

C = kat\\y\\p 

D = (1 - (1 - AA l3 )~ 1 CA l3 )~ 1 
F = DB(1 -AA^A?. 



\\x\\ s , t ,oo [l-kM^l-AA^^A?] 
< |x s |+A : ao||y||Ml-^A' 3 )- 1 A /3 . 



sup |a? r | < (1 - (1 - AA' 3 )- 1 CA /3 )~ 1 

0<r<t 



sup \x r \ +B(1 -AA /3 )- 1 A /3 

0<r<s 



Or 



Denote 



sup lav | < D sup \x r \ + F. 

0<r<t 0<r<s 



Z n = sup |x r |, 

0<KnA 



where n = ? . Then 



Z n < DZ n _x + F < ■ ■ ■ < D n Z + DkF - 



k=0 
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This yields 



sup |a; t | < (l-il-AA^CA^-^xol 

0<t<T 

n-1 



+ 2(1 - (1 - AA /3 )- 1 CA /3 )- fe - 1 B(l - AA^A 13 . 



fe=0 



Then we let A satisfy 



AA 13 < 1/3, CA P < 1/3, BA 13 < 1/3 



Namely, we take 



1 



3(AV5VC) 



1//3 



Then 



sup H < 2 T / A (|a; | + 1) 

0<t<T 



< 2 feT[||/'|UVa Va 1 ] 1/,3 ||y|l^(| Xo | + !) . 



The proof of the theorem is now complete. ■ 

Suppose now that we have two differential equations of the form 



x t =x + / f(x s )dy s 
Jo 



x t = x + / f(x a )y a , 



and 



where y and y are Holder continuous functions of order (3 > 1/2, and / and / 
are two functions which are continuously differentiable with Holder continuous 
derivatives of order A > jj — 1. Then, we have the following estimate. 

Theorem 3.2 Suppose in addition that f is twice continuously differentiable 
and f" is bounded. Then there is a constant k such that 

Sup \x r -X r \ < P Wvt^T 

0<r<T 



X< \x - X \ + \\y\\o,T,0 



||/-/||oo + |k||0,T, /3 ||/ / -/ / | 



+ H/lloo + ||/||oolMlo,r,oo \\y-y\\o,T,0 



where 



D = H/'IU V (II/'IIooIIi/IIct,/? + ||/"||oo(||x||o^ + \\x\\o,t, )T^ 
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Proof. Fix s,t e [0,T]. Set 

x t -x t - (x s - x s ) = h+ I 2 + h 

where 



h = I [f{x r ) - f(x r )]dy r 

J s 



h = 
h = 



dy r 



f(x r ) - f(x r ) 

/ f{x r )d[y r -y r ]. 

J S 

The terms I2 and I3 can be estimated easily. 

|/ 2 | < k\\ y y [\\f f\ut sf + - fuwxUtAt - ^ 

and 

|/ 3 | < k\\y vh [ll/lloo(i - sf + \\~f\U\xUAt - s) W , 

where \\y\\p = \\y\\o,T,p and \\y-y\\p = \\y-y\\o,T,p- The term h is a little more 
complicated. 

I'll < f s \D^ + [f{x r )-f{x r )]\\D]z a y t -{r)\dr 

< k f hUAt - r)^- 1 [\f(x r ) f(x r )\(r s)- a 

J s 

+ \\}'\\oc\\x - Z\\ 8 ,rAr - sY"* 

+ ||/"||oo|k - £|| a , r ,oc [\\x\\s,r,P + \\x\\,,r,p] (r - S^""] dr 

< *||l/||/J {H/'IUI* ~ 5|| a ,t,oo(* - S)' 3 + H/'Hoolk - i|| a , t>/J (t - S ) 2/3 

+ Iiri|oo||.x - i||,, t>0 c [||x|| M/J + \\x\\ B , ttP ] (t - s) 2f3 } . 

Therefore 

\\X - X\\s,t,0 < k\\y\\p {||/'||oo|k - 5||a, t ,oo + ||/'||oo||z - x\\ s , tt0 (t - sf 

+ ||/"|| oc ||x - jc|| aiti00 [\\x\\ a , tt p + \\x\\,,t,fi] (t - sf 

+ ||/-/||oo + ||/'-/'||ooP||, > t,/3(i-^} 

+k\\y-n0[\\f\\oo + \\f'\\oo\\x\\s,tAt-sf] ■ 
Rearrange it to obtain 



\\x-x\Ui3 < m-kinuiyWf^t-sf)- 1 ! \\ y y 
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+||/"||oo||a; - x\\s,t,oo [\\x\\s,t,p + \\x\\s,t,(}] (t - s) 13 

+ \\f-f\\oc + \\f-f'\U\x\\s,tAt-s) 
+fc||l/-»ll/j[ll/l|oo + ||/'||oo||i|| a ,t,|8(*-*) 

Set A = t - s, and A = fc||/'||oo||y||/3- Then 

\\x - x\\s,t,oo < \x s - x 3 \ + \\x - x\\ a ,t,p{t - s) 
< \x,-x,\+k(l- AA /3 )- 1 A /3 

+ ||/"|| oc ||x - x\\ att>00 [\\x\\ a>t ,p + \\x\\,,t,p] A^ 

+II/-/Hoo + ||/'-/'||ooPIIm,/3A' 3 



% %\\s,t,c 



+k\\y-yh 



oo + 11/ ||oo|p|| s ,t,/3 



A' 3 



Denote 



Then 



B = k\\y\y (ll/'IU + ||/"||oc(||a:||o,r,/ 3 + \\i\\o,T,p)T^) . 

iiz-xii,, t>0 o < (l-a-^r^B) -1 

x{|x s -i s | + fc(l-AA' 3 )- 1 A' 3 

||2/||J||/-/||oo + ||/'-/'||ooP|| S> t l/3 A^ 



+ 511/8 



A" 



Let A satisfy 
Namely, we take 

Then 



AA^<l/3, BA fi <l/3 



1 

3 (A V S) J 



I a; - x\\ s ,t,oc < 2 [|z s -x s \ + CA ] 
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where 



C 



Wvh ll/-/llcx 3 + ]]/'-/ / IUI|x|| s ,^ 



+\\v-vh 



ll/IU + ll/'lloollilU,*,^' 



Applying the above estimate recursively we obtain 

sup \x r - x r \ < 2™ [|xq - x \+ CA' 3 ] , 

0<r<T 

where T = nA. Or we have 

sup \x r -x r \ < 2 fe(ll/ ' IUv ( |l/ ' lt ~ + ll/ '' IU(l|x|lo ' T ^ +l|5|lo ' T ^ )T3 ) )1/ ^ l|y|| o^^ T 

0<r<T 



x\\x - xa\ + \\y\\o,T,p 11/ - /Hoc + \\x\\o,T,0\\f - fl 



+ ll/lloo ||£c||o,t,oo \\y - y\\o,T,p 



4 Stochastic differential equations driven by a 
fBm 

Let B — {B t ,t > 0} be an m-dimensional fractional Brownian motion (fBm) 
with Hurst parameter H > 1/2. That is, B is a Gaussian centered process with 
the covariance function E{B l t B{) = Rii(t,s)Sij, where 

R H (t,s) = ~(t 2H + s 2H -\t-s\ 2H ). 
Consider the stochastic differential equation 



X t — x 



a(X s )dB s 



(4.1) 



This equation has a unique solution (see and [Hj) provided a is continuously 
differentiable, and a' is bounded and Holder continuous of order A > -4 — 1. 
The stochastic integral is interpreted as a path-wise Riemann-Stieltjes integral. 

Then, using the estimate ()3.4(l in Theorem 13. II we obtain the following esti- 
mate for the solution of Equation (S3) , if we choose (3 £ Notice that 



A < 2 
< - 



sup |X t | < 2 fcT (ll CT 'll~ v l^l)ll<^(|X | + 1) , 

0<t<T 



(4.2) 
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If <j is bounded we can make use of the estimate (I3.2|) and we obtain 



sup \X t \ < |X | + fcT||a|| oo ||<T'||^|| J B||J T/3 . (4.3) 

0<t<T ' 

These estimates improve those obtained by Nualart and Ra§canu in j3] based 
on a suitable version of GronwalPs lemma. The estimates (|4.2|) and (|4.3|l allow 
us to establish the following integrability properties for the solution of Equation 

Theorem 4.1 Consider the stochastic differential equation \4- If a' is bounded 
and Holder continuous of order A > -h — 1, then 

E ( sup \X t \ p ) < oo (4.4) 

\0<t<T J 

for allp>2. If furthermore a is bounded, then 



.E^expA^sup \X t pjj <oo (4.5) 

for any A > and 7 < 2/3. 

If we apply these results to the linear equation satisfied by the derivative in 
the sense of Malliavin calculus of X t then we get that X t belongs to the Sobolev 
space IB> 1 ' P for all p > 2. This implies that if the coefficient a is infinitely 
differentiable with bounded derivatives of all orders, then, X t belongs to D°°. 
This allows us to deduce the regularity of the density of the random vector X t 
at a fixed time t > assuming the following nondegeneracy condition: 

(H) The vector space spanned by j (cr lJ {X Q )) 1<l<d , 1 < j < m\ is R m . 

That is, we have the following result. 



Theorem 4.2 Consider the stochastic differential equation \4-l\l - Suppose that 
a is infinitely differentiable with bounded derivatives of all orders, and the as- 
sumption (H) holds. Then, for any t > the probability law of X t has an C°° 
density. 

In 0j Nualart and Saussereau have proved that the random variable X t 
belongs locally to the space HJ°°, and, as a consequence, they have derived the 
absolute continuity of the law of X t under the assumption (H) . 
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